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Abstract: We introduce a supersymmetric Chern-Simons theory whose low energy 
physics is that of the fractional quantum Hall effect. The supersymmetry allows us to 
solve the theory analytically. We quantise the vortices and, by relating their dynamics 
to a matrix model, show that their ground state wavefunction is in the same universality 
class as the Laughlin state. We further construct coherent state representations of 
the excitations of a finite number of vortices. These are quasi-holes. By an explicit 
computation of the Berry phase, without resorting to a plasma analogy, we show that 
these excitations have fractional charge and spin. 
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1. Introduction and Summary 

Supersymmetry is a much beloved tool of high energy theorists. Supersymmetric field 
theories are often tractable, even at strong coupling, yet remain rich enough to exhibit 
a wide range of interesting dynamics. 
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In contrast, supersymmetric theories are much less studied in the condensed matter 
community, even in the limited role of toy models for strongly coupled phenomena. In 
part this is because supersymmetry typically provides analytic control for relativistic 
theories at vanishing chemical potential. At finite density, where most problems of 
interest in condensed matter lie, supersymmetry is usually broken and any advantage 
it brings is lost 1 . 

There is, however, a class of theories in d — 2 + 1 dimensions which are supersym¬ 
metric, yet non-relativistic [7]. In these theories, supersymmetry is retained even at 
finite density. Despite the vast literature on supersymmetric field theories, the quantum 
dynamics of these models remains relatively unexplored. The purpose of this paper is 
to show that the low-energy physics of these theories is that of the fractional quantum 
Hall effect. 

Of course, the fractional quantum Hall effect is one of the most studied topics in 
physics over the past three decades. The theory rests on a beautiful and intricate web 
of ideas involving microscopic wavefunctions [8], low-energy effective Chern-Simons the¬ 
ories [9, 10, 11, 12, 13] and boundary conformal theories [14, 15]. The supersymmetric 
theory that we present here is unlikely to be of direct relevance to any material. In¬ 
stead, it should be viewed as a toy model whose role is to highlight some of the links 
between these different approaches to the quantum Hall effect. 

In the rest of this introduction, we describe the supersymmetric model in more 
detail and explain what it’s good for. It is an Abelian Chern-Simons theory, coupled 
to both bosonic and fermionic non-relativistic matter fields. In this manner, it is an 
amalgamation of effective theories of [11] and [12]. The model has vortices and these 
are viewed as the “electrons”. The vortices are “BPS objects” [16]: this means that 
they experience no classical static forces. It also means that they are protected by 
supersymmetry in a way which we describe in the main text. This property allows us 
to perform an explicit quantisation of the vortex dynamics. We show that the ground 
state wavefunction of the vortices lies in the same universality class as the Laughlin 
wavefunction. It has the same long range correlations, but differs on short distance 
scales. 

We also describe the excitations of a droplet of vortices. There are gapless, chiral 
edge excitations which, we show, are governed by the usual action for a chiral boson 

Where are a number of notable exceptions, including the role of supersymmetry in disorder [1], the 
possibility of emergent supersymmetry [2, 3, 4, 5] and the study of supersymmetry protected phases 
[ 6 ], 
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[17], suitably truncated due to the presence of a finite number of vortices. Finally, 
we construct the quasi-hole excitations in this model and compute their Berry phase. 
This is, of course, a famous computation for the Laughlin wavefunctions [18]. However 
the usual analysis relies on the plasma analogy [8], and the (admittedly well justified) 
assumption that the classical 2d plasma exhibits a screening phase. In contrast, here 
we are able to perform the relevant overlap integrals analytically to show that the 
quasi-holes have the expected fractional charge and statistics. 

Many of the properties of vortices described above follow from the fact that their 
dynamics is governed by a quantum mechanical matrix model, which was introduced by 
Polychronakos to describe quantum Hall physics [19] and further studied in a number 
of works [20, 21, 22, 23]. The results of this paper show how this matrix model is 
related to more familiar effective field theories of the quantum Hall effect. 

The paper is organised as follows. In Section 2 we introduce the non-relativistic, 
supersymmetric theory. After a fairly detailed description of the symmetries of the 
theory, we discuss its two different phases and its spectrum of excitations. Section 3 is 
devoted to a study of BPS vortices and contains the meat of the paper. We will show 
that the low-energy dynamics of vortices is governed by the matrix model introduced 
in [19]. We review a number of results about this matrix model and derive some new 
ones. Finally, in Section 4 we offer some ideas for the future. A number of calculations 
are relegated to appendices. 

2. Non-Relativistic Chern-Simons-Matter Theories 

We start by introducing the d = 2 + 1 non-relativistic, supersymmetric Chern-Simons 
theory. The theory consists of an Abelian gauge field A^, coupled to complex scalar 
held 0 and a complex fermion i/j. The action is 

S = I dtd 2 x !\i(p ] D o0 + i^Vo'ip - V «' , P ~ ~^ e ^ P A ^» A p 

-Mo + ] Bif) - (|0| 4 - Ml 2 + 3 I0| 2 |'0| 2 ) ) (2- 1 ) 

2 m mk J 

Some conventions: the subscripts p — 0,1, 2 run over both space and time indices, 
while a = 1,2 runs over spatial indices only. The fermion carries no spinor index. Both 
4> and ip are assigned charge 1, so the covariant derivatives read = d^ — iA^cj) and 
similarly for %b. The magnetic field is B — d\ A 2 — c^Tli- Finally j^l 2 = VM — — 
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There are three parameters in the Lagrangian: the Chern-Simons level k G Z + , the 
mass m of both bosons and fermions, and the chemical potential /i. As we will see later, 
the chemical potential /j can be more fruitfully thought of as a background magnetic 
field for vortices. 

The first order kinetic terms mean that the action (2.1) describes bosonic and 
fermionic particles, but no anti-particles. The quartic potential terms correspond to 
delta function contact interactions between these particles. In the condensed matter 
context, the gauge field is considered to be emergent. One of its roles is to attach flux 
to particles through the Gauss’ law constraint, which arises as the equation of motion 
for A 0 , 


b= T w 2 +m 2 - q (2.2) 

We’ll learn more about the importance of this relation later. 

The action (2.1) can be constructed by starting from a relativistic Chern-Simons 
theory with J\f = 2 supersymmetry and taking a limit in which the anti-particles 
decouple. For the case n = 0, this was first done in [7] and we review the procedure in 
Appendix A. To our knowledge, the supersymmetric theory with /i / 0 has not been 
previously constructed, although the bosonic sector of our theory is similar, but not 
identical, to a model studied by Manton [24] which shares the same vortices as (2.1). 
We will describe these vortices in some detail in Section 3. 

2.1 Symmetries 

The action (2.1) is invariant under several symmetries. A number of these play an 
important role in what follows and we provide the relevant details here. 

Bosonic Symmetries 

Invariance under time translations gives rise to the Hamiltonian. After imposing the 
Gauss’ law constraint (2.2), this takes the concise form 

H= - [ d 2 x \V z <t>\ 2 + \V- Z ^\ 2 + j\(t>\ 2 \^\ 2 (2.3) 

m J k 

where z = x 1 + ix 2 and z = a; 1 — ix 2 . Correspondingly, d z = |(<9i — id 2 ) and d z = 
\{d\ + id 2 ). 
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Invariance under spatial translations gives rise to the complex momentum, P = 
\{P\ — 1 P 2 ), which we write as 


P = 



d 2 x zB with P 


d 2 x ^V z <j> - V 


(2.4) 


The P contribution is the standard Noether charge for spatial translations. The second 
term, proportional to the chemical potential /i, requires some explanation. As shown in 
[25], it arises because a translation is necessarily accompanied by a gauge transforma¬ 
tion so that, for example, 5$ = The presence of the chemical potential term /jAq 
in the action then means that naive Noether charge for translations is not gauge invari¬ 
ant. This is remedied by the addition of a total derivative, resulting in the improved, 
gauge invariant momentum above. Note, however, that the resulting momentum P is 
not itself translationally invariant. We shall comment further on this below. 


A similar subtlety occurs for rotations. The conserved angular momentum is given 
by 


J = I d 2 x (z$V z (t) - zV z tf 0 + zip'V z ip - zV- z 0 - - |N 2 ^) (2.5) 

Here the first terms are standard. The penultimate term, ]0| 2 , arises because the 
fermion is taken to have spin 1/2. The final term again arises as an improvement term 
in the Noether procedure which ensures that the resulting angular momentum is gauge 
invariant [25]. 

The number of bosons and fermions in our model are individually conserved. The 
corresponding Noether charges are 

Mb = J d 2 x (^(j) and Mf = J d 2 x (2.6) 

The total particle number is simply the charge under the Abelian gauge group 

M = Mb Mf 

We denote the axial combination as 

TZ = Mb — Mf 

This will play the role of an R-symmetry in the supersymmetry algebra. 
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The presence of the anomalous term in the expression for the momentum (2.4) has 
an interesting effect on the commutation relations. (Here we describe the quantum 
commutation relations rather than classical Poisson brackets.) We find 

[H,P] = ~^P and [H,P]= 0 (2.7) 

TflKi 

So the Noether charge P is conserved, but the translationally invariant momenta P 1 
and P act as raising and lowering operators for the spectrum. Further, the conserved 
momenta do not commute. We have 

{P,P<]=-^M (2.8) 

Both (2.7) and (2.8) are similar to the commutation relations in quantum mechanics 
for momenta in a magnetic held. This is because, as we will describe in more detail 
below, /i acts like an effective magnetic for vortices while the Gauss’ law constraint 
ensures that all excitations carry some vortex charge. 

When fi = 0 the theory also enjoys both a Galilean boost and, more surprisingly, a 
(super)conformal symmetry [7] (see also [26]). In this paper, we restrict ourselves to 
the non-conformal theory with chemical potential fi > 0. 

Supersymmetries 

The action (2.1) enjoys two complex supersymmetries [7]. These often go by the name 
of kinematical and dynamical supersymmetries. The kinematical supersymmetry is the 
simpler of the two. 

5i(f) = e\ip , <5i ij) = -ei (p , S±A Z = 0 , 5\A 0 = (ei<ppA - e\i>(p ] 

This is reminiscent of the Green-Schwarz spacetime supersymmetry on the string world- 
sheet. The transformation on (p and ip is generated by 

Qi — J d 2 x (pUp (2.9) 

This does not specify a transformation for Ho which does no harm as long as we allow 
ourselves to impose Gauss’ law. We will see the implications of this below. 

Under the dynamical supersymmetry, the fields transform as 

d 2 (p = e\V z ip , 8 2 ip = e 2 V z (f> , 5 2 A Z = -^el'ippA 
8 2 Aq = (U 2 (fUv z tp - e 2 (pT> z ip 
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with supercharge 


Q 2 = j d 2 x <pV- z i\> (2.10) 

The supersymmetry algebra is 

{QuQ\} = M , IQ2,Q i 2 } = jH , {Qi,Ql}=P (2.ii) 

Note that the two supercharges generate the translationally invariant momentum P, 
rather than the conserved momentum P. There is also a mild surprise in the commu¬ 
tators of bosonic and fermionic charges 2 , in particular 

\n, Q I } = ~Q l ( 2 . 12 ) 

This means that although the kinematic supersymmetries leave the action invariant, 
when /i / 0 they do not result in a symmetry of the spectrum. This can be traced to 
the fact that Gauss’ law was required, both in the construction of the Hamiltonian (2.3) 
and in the derivation of the commutators (2.12). (To our knowledge, the possibility of 
non-relativistic supersymmetry generator Q i was first raised in [27] who also pointed 
out that this generator is spontaneously broken in any vacuum with non-vanishing 
particle number.) Other commutators follow from Jacobi identities and give [Q 2 ,H] = 
[Qi,P] = [Q\,P] = 0 while [ Q 2 ,P] = (2/m)[H,Q 1 \. 

Finally, the commutators of the angular momentum will also be important for our 
story. There’s nothing unusual about them. We have 

[J,Qi] = ~Qi and [J,Q 2 ] = ^Q 2 (2.13) 

which is the statement that the supercharges have spin =pl/2. The total particle number 
commutes with all supersymmetries, [AT , Q a ] = 0, but only because of cancellations 
between boson and fermion numbers. Individually, we have 

Qa] Qa and [Afp; Qa\ T Qa 

This justifies our previous claim that TZ is an “R-symmetry” since [JZ, Q a \ = 2 Q a . 
From these, we deduce that 

[J-\Mf,Q2} = 0 (2.14) 

This fact will be important in Section 3.2. 

2 We thank Nima Doroud for very useful discussions regarding this algebra. 
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2.2 The Vacuum, The Hall Phase, and Excitations 

We now describe some basic features of the dynamics of our model. Because non- 
relativistic field theories have no anti-particles, the theory decomposes into sectors 
labelled by the conserved particle numbers which, in our case, are Mb and Mf- To 
solve the theory, we need to determine the energy spectrum in each of these sectors. 

One way to organise these sectors is to start with the M = 0 Hilbert space and build 
up by adding successive particles. Instead, we will take a dual perspective. Our theory 
enjoys a topological current, 

J' 1 = ^-e lwp d,A P (2.15) 

Z IT 

The associated particles are vortices. We will view these vortices as the “electrons” of 
our theory. 

Our theory has two translationally invariant ground states consistent with Gauss’ 
law (2.2), both of which have H = 0. We call these the vacuum and the Hall Phase. 
They are given by 

The Vacuum: |0| J = /j and B — 0 (2-16) 

The Hall Phase: \<b \ 2 = 0 and B =- -J— (2.17) 

k 

The vacuum state contains no vortices, f d 2 x J° = 0. However, the bosons have 
condensed which means that the particle number is M = oo. In contrast, the Hall 
phase has vanishing particle number but infinite vortex number, J d 2 x J° — oo. 

The purpose of this paper is to understand what happens as we inject vortices into the 
vacuum. For any finite number of vortices, the system breaks translational invariance. 
But, as we fill the plane with vortices, the Hall phase emerges. In Section 3, we tell 
both the classical and quantum versions of this story in some detail. First, however, 
we describe some simple properties of excitations above each of these ground states. 

The Vacuum 

The key feature of the vacuum state is that f/(l) gauge symmetry is broken. This 
ensures that the theory admits topological, localised vortex solutions. These vortices 
will be the main focus of this paper and we postpone a more detailed discussion until 
Section 3. Here, we summarise their three main properties: 

• Vortices are gapless. States with an arbitrary number of vortices exist with H = 0. 



• Vortices have statistical phase nk. This means that the vortices are bosons when 
k is even and fermions when k is odd. 

• Vortices are singlets under supersymmetry. 

There are further excitations above the vacuum arising from the fundamental holds 
0 and 0. These excitations are both gapped, with an excitation energy 27Tfi/mk. These 
excitations can be generated from the vacuum by using the raising operators P t and 
Q\, together with the supercharges Q 2 and Ql- 

The Hall Phase 

The Hall phase has an unbroken [/(1) gauge symmetry and the long-distance physics is 
dominated by the Chern-Simons term, ft is well known that such theories capture the 
essential properties of the fractional quantum Hall effect. We now take the opportunity 
to review this standard material (see, for example, [28, 29] for reviews). 

To describe quantum Hall physics, it is not enough to specify the Lagrangian; we 
need to know how electromagnetism couples to the theory 3 . (Recall that the Abelian 
gauge held A p in the Lagrangian (2.1) should be thought of as an emergent, statistical 
gauge held, not the electromagnetic held.) Since we wish to treat the vortices as the 
“electrons” of the theory, the background electromagnetic held A® xt must couple to the 
topological current (2.15), 

£ Ha ii = ^e^A,d u A p + eAfJ* + ... 

Here e denotes the electron charge, while ... includes the rest of the Lagrangian (2.1), 
as well as the d — 3 + 1 dimensional Maxwell term for kL xt . 

A 1 

We momentarily ignore the fundamental helds 0 and %b. Integrating out A p , the 
quadratic Lagrangian for the background held is given by 

£ HaH + - 

The effective action S' e fr[A ext ] = J d 3 x Cna\h is now a functional of the non-dynamical, 
background electromagnetic held. Its role is to tell us how the system responds to an 

Where are two, dual, descriptions of the long-wavelengtlr quantum Hall physics in terms of Chern- 
Simons theories. In one description, the Chern-Simons level is equal to is, the filling fraction [9, 10], 
the electrons are the fundamental excitations and the vortices the fractionally charged quasi-particles. 
Here we are interested in the dual description, related by a particle-vortex duality transformation, 
where the Chern-Simons coefficient is 1/v and the electrons are vortices. 
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applied electromagnetic field through the relation (J M ) = dS e g /cM® xt . The result is a 
Hall conductivity, 


a H = 


2nk 


(2.18) 


This is the response of a fractional quantum Hall fluid at filling fraction v = 1/k. 


Let us now return to the fundamental fields 0 
and 0. Each of these experiences a magnetic held 
B = —27 T/x/k and forms Landau levels. The usual 
Landau level quantisation results in a spectrum 

Ell = -—~(l + 1 / 2 ) 

m 

with l = 0,1,.... However, the Lagrangian (2.1) 
also includes extra terms which shift the over¬ 
all energy of these states. The shift is down for 
bosons and up for fermions, as shown in the fig¬ 
ure. The net result is that the energies of the 
Landau levels, at leading order, are given by 


Energy 


Bosons 

/ 

\ 

Fermions 







"4 




-3 







-2 







-1 








Figure 1: Fundamental excitations 


2tc/lI I l — 0,1, 2,... for 0 
1 l — 1, 2,... for 0 


The gapped states, with l > 1, arising from 0 have spin 1/2 k\ those arising from 0 
have spin (1 + k)/2k. The Gauss law constraint (2.2) ensures that, when coupled to 
a background electromagnetic held, each of these carries charge —e/k. These are the 
quasi-particle excitations of our supersymmetric quantum Hall fluid. The supercharges 
Q 2 and Q\ map between the fermionic and bosonic gapped Landau levels. 

The system also has an a gapless band of quasi-particles, arising from the lowest 
Landau level of 0. These modes are not free: they interact through the 0 4 potential 
in (2.1). Nonetheless, supersymmetry ensures that these states have vanishing energy 
at all orders of perturbation theory. This is because the commutation relations for Q 2 
require that any excitation with H > 0 must be paired with an excitation that differs 
by spin 1/2. Yet the states in lowest Landau level have no partners and must, therefore, 
remain at zero energy. In essence, the theory has an infinite Witten index Tr(—1) F . If 
we start from the lowest Landau level, we can build up to higher levels by acting with 
Pt and Q\. 
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The presence of a gapless Landau level may appear to contradict our claim that 
this system describes quantum Hall physics. After all, one of the defining features of 
a quantum Hall state is that it is gapped and incompressible. We will resolve this 
in Section 3 by studying how the Hall phase emerges from vortices when placed in a 
confining potential. We will show that, for any finite number of vortices, there is a 
unique incompressible droplet of lowest angular momentum. However, in the absence 
of a confining potential, this droplet has zero energy edge modes and zero energy 
quasi-hole excitations. The gapless Landau level describes these degrees of freedom for 
an infinite number of BPS vortices, an interpretation recently suggested in a different 
context in [30]. We will revisit this in Section 3.4 in the context of the non-commutative 
approach to quantum Hall physics. 

ft is worth mentioning that this situation is not unusual in quantum Hall systems. 
The special, ultra-local Hamiltonians (such as Haldane pseudo-potentials) used as mod¬ 
els of quantum Hall physics also have zero energy edge modes and zero-energy quasi-hole 
excitations for finite droplets. See, for example, [31, 32] for related discussions. 


3. A Quantum Hall Fluid of Vortices 


We would like to understand how to interpolate from the vacuum to the Hall phase. 
We do this by injecting vortices. These vortices are BPS which, in this context, means 
that they have H = 0 and lie in a protected sector of the theory. From the form of the 
Hamiltonian (2.3) and Gauss’ law (2.2), it is clear that solutions with vanishing energy, 
H — 0, can be constructed by solving the equations 

V z 4> = 0 and B = — (|0| 2 — n) (3.1) 

k 

with the fermions set to zero: = 0. 


The vortex equations (3.1) are well studied. Solutions are labelled by the integer 
winding of the scalar field (j) or, equivalently, by the magnetic flux 


1 


n — — 


2vr 


d 2 x B e Z+ 


(3.2) 


In the sector with winding n, the most general solution to (3.1) has 2 n real parameters 
[33, 34], These parameters are referred to as collective coordinates or, in the string 
theory literature, moduli. When vortices are well separated, these correspond to n 
positions on the complex plane. The existence of these moduli reflects the fact that 
the coefficient of the quartic interaction in (2.1) has been tuned to the critical value, 
ensuring that there are neither attractive nor repulsive forces between the vortices. 
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Figure 2: Two points in the moduli space of n = 7 vortices 


As vortices coalesce, they lose their individual identities and the interpretation of 
these moduli changes. It is tempting to label the vortex by the point at which the 
Higgs field vanishes, but this does not provide an accurate description of what the 
vortex profile looks like. Instead, as we show in Section 3.4, in this regime it is better 
to think of the 2 n moduli as describing the edge modes of a large, incompressible fluid. 

Why do Vortices Form a Fractional Quantum Hall State? 

The rest of this section is devoted to a detailed analysis of the quantum dynamics of 
vortices. We will ultimately show that their ground state is given by the Laughlin 
wavenfunction. But here we first provide a hand-waving argument for why we expect 
the vortices to form a quantum Hall fluid. 

We first note that the chemical potential term (j,Aq, present in the Lagrangian (2.1), 
can be viewed as a background magnetic field for vortices. It can be written as 



where J M is the topological current (2.15) and 


nex t 

^ =-With 


^ext = ^ 


(3.3) 


e 


This means that we expect the dynamics of vortices to correspond to particles moving 
in a background magnetic field. Nonetheless, it may be rather surprising that the 
vortices form a Hall state because, as we have seen, there is no force between the 
vortices. Yet the key physics underlying the fractional quantum Hall effect is the 
repulsive interactions between electrons, opening up a gap in the partially filled Landau 
level. 

Although there is no force between vortices, they are not point particles. Instead, 
they are solitons obeying non-linear equations and, as they approach, the solutions 
deform. Indeed, when the vortices are as closely packed as they can be, they form 
a classically incompressible fluid as shown in the right-hand side of Figure 2. The 
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scalar field 0 has an n th order zero in the centre of the disc and numerical studies show 
that the solution is well approximated as a disc of magnetic flux in which 0 = 0 and 
B = —27 rp/k. This has motivated the “bag model” of vortices in [35, 36]. For us, it 
means that the vortex is a droplet of what we have called the “Hall phase”. 


When n vortices coalesce, the radius R of the resulting droplet can be estimated 
using the flux quantisation (3.2) to be 


R 



(3.4) 


Now we can do a back-of-the-envelope calculation. In a magnetic field B ext , the number 
of states per unit area in the lowest Landau level is eB ext /2ir = p. In an area A = 
nR 2 = nk/p , the lowest Landau level therefore admits B ext A = nk states. We’ve 
placed n vortices in this region, so the filling fraction is 


v = 


1 

k 


This, of course, is the expected filling fraction in the Hall phase with conductivity 
(2.18). 


eB ext 

2vr 


b 


3.1 The Dynamics of Vortices 


7 tR 2 = 


nk 

b 


We now turn to a more detailed description of the dynamics of vortices. We first intro¬ 
duce the vortex moduli space, A4 n . This is space of solutions to the vortex equations 
(3.1) with winding number n. As we have already mentioned, 


dim(Al n ) = 2 n 


The coordinates X a , a = 1,... ,2 n, parameterising A4 n are the collective coordinates 
of vortex solutions: (p(x] X) and A a (x\X). 

The standard approach to soliton dynamics is to assume that, at low-energies, motion 
can be modelled by restricting to the moduli space [37]. This is usually applied in rela¬ 
tivistic theories where the action is second order in time derivatives and typically pro¬ 
vides an accurate approximation to the real dynamics. Here we have a non-relativistic 
theory, first order in time derivatives, and this results in a number of differences which 
we now explain. One ultimate surprise — which we will get to in Section 3.2 — is 
that there is no approximation involved in the moduli space dynamics in this system; 
instead it is exact. 
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The first, and most important difference, is associated to the meaning of the space 
A4 n . In relativistic theories, JA n is the configuration space of vortices and the dynamics 
is captured by geodesic motion on A4 n with respect to a metric g ab (X). It is known 
that A4 n is a complex manifold, with complex structure J, and the metric g a b(X) is 
Kahler. For completeness, we explain how to construct this metric in Appendix B. 

In our non-relativistic context, it is no longer true that A4 n is the configuration space 
of vortices. Instead, it is the phase space. The dynamics of the vortices is described by 
a quantum mechanics action of the form, 

Cortex = Jdt X a (X)X a (3.5) 

where X{X) is a one-form over A4 n . Our goal is determine this one-form. 

In fact, this problem has already been solved in the literature. A model which shares 
its vortex dynamics with ours was previously studied by Manton [24] and subsequently, 
in more geometric form, in [38, 39]. The main result of these papers is that T is an 
object known as the symplectic potential. It has the property that 

dT = 12 (3.6) 

where 11 is the Kahler form on A4 n , compatible with the metric g a b and the complex 
structure J. 

For a single vortex, the moduli space is simply the plane C and the Kahler form is 

12 = — dz A dz 
2 

For n > 2 vortices the Kahler form is more complicated. We describe the construction 
of 12 in Appendix B. Explicit expressions are only known for well-separated vortices 
[39]. 

The derivation of (3.6) given in [24, 38, 39] relics on a parameterisation of the vortex 
moduli space introduced earlier in [41]. The use of these coordinates means that the 
calculation is not entirely straightforward. For this reason, in Appendix B, we present 
a simpler derivation of (3.6) which does not rely on any choice of coordinates. (For a 
different approach to particle dynamics appropriate for vortices, see [40].) 
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There are no Fermion Zero Modes 

The vortices are BPS states: they are annihilated by the supercharge Q 2 . In the context 
of first order dynamics, this means that the collective coordinates X do not transform 
under Q 2 . In particular, there are no accompanying Grassmann collective coordinates. 
Indeed, it is simple to check explicitly that there are no fermionic zero modes in the 
background of the vortex. 

The upshot is that the vortices themselves are supersymmetric singlets. The role of 
supersymmetry in the vortex dynamics (3.5) is restricted to ensuring the vortices have 
strictly vanishing energy, H = 0, in the full quantum theory. 

The fact that the BPS solitons have no fermion zero modes may come as something 
of a surprise. Indeed, it is rather different from what happens for BPS solitons in rela¬ 
tivistic field theories or in string theory. It is worth pausing to explain this difference. 
In more familiar relativistic theories, if a soliton is invariant under a given supercharge 
Q then that supercharge will descend to the worldvolume theory, relating bosonic and 
fermionic zero modes on the worldvolume. However, when we say that a soliton is 
invariant under Q, we mean that the static configuration is invariant: when the soliton 
moves, the supercharge Q typically acts and generates a fermionic zero mode. This 
means that while Q does not act on the bosonic configuration space of the soliton, it 
does act on the phase space. 

In our non-relativistic theory, the statement that Q 2 annihilates the soliton is stronger: 
it means that Q 2 does not act on the soliton phase space. This is the reason that there 
are no associated fermionic zero modes. 

3.2 Introducing a Harmonic Trap 

We have derived a low-energy effective action (3.5) for the vortex dynamics. However, 
this dynamics is boring. The equation of motion arising from (3.5) is 

n ab x b = 0 => X a = 0 

The vortices don’t move. They are pinned in place. 

The lack of dynamics follows because there is no force between vortices and, in a 
first order system, we don’t have the luxury of giving the vortices an initial velocity. 
To get something more interesting, we impose an external force on the vortices. We 
will do so by introducing a harmonic trap. We want this trap to be compatible with 
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supersymmetry. We can do this by choosing the new Hamiltonian 

^new = H + t) ^ J ~ ~MlF 

where J is the angular momentum (2.5), Mf the fermion number operator (2.6) and 
to dictates the strength of the trap. From (2.14), we see that this Hamiltonian remains 
invariant under Q 2 , although not Q\. When evaluated on BPS vortices, the Hamiltonian 
is simply 

tfnew = ~ J dM \Z\ 2 B (3.7) 

This new Hamiltonian is the angular momentum of a given BPS vortex configuration: 
it preserves the BPS nature of vortices while shifting their energy. Evaluating (3.7) 
on a vortex configuration provides a function J(X) over the vortex moduli space A4 n 
which governs the their low-energy dynamics, 

Svortex = f dt ( 'r a (X)X a -uJ(X )) (3.8) 

We will now look at some examples of the classical dynamics dcescribed by this action. 

Classical Motion in the Trap 

The harmonic trap (3.7) favours those vortex solutions that are clustered towards the 
origin. The lowest energy configuration now has all vortices coincident at the origin, as 
in the right-hand picture in Figure 2. As we have seen, the size of this vortex is given 
by (3.4), so the angular momentum of this state is 

Jo ~ —| / dr 2nr 3 B = (3.9) 

This is the only static configuration. All other solutions evolve through the equation 
of motion 

(9 7 

n ab x b = tu-xM (3.10) 

In particular, a single vortex displaced a distance r 3> l//i from the origin, will have 
angular momentum J fa 7r/xr 2 . This vortex orbits around the origin with frequency to. 

There is something rather surprising about the moduli space approximation for this 
first order dynamics: it is exact! The solutions to the equation of motion in the presence 
of the trap are simply time dependent rotations of the static solutions so, for example, 
(f) = <f>(x;X(t)), with X(t) obeying (3.10). This a property of any first order system 
with a Hamiltonian, such as H = J, which acts as a symmetry generator on the moduli 
space. 
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3.3 The Quantum Hall Matrix Model 

The description of the vortex dynamics (3.8) is, unfortunately, rather abstract. For 
n > 2 vortices, we have only implicit definitions of the Kahler form fl and the angular 
momentum J on the vortex moduli space. It seems plausible that one could make 
progress using the parameterisation of the vortex moduli space introduced in [41]. 
Here, however, we take a different approach. 

An alternative description of the vortex moduli space is provided by D-branes in 
string theory [42], This is analogous to the ADHM construction of the instanton moduli 
space. The vortex moduli space M n is parameterised by: 

• An n x n complex matrix Z 

• A n-cornponent complex vector p 

These provide n(n + 1) complex degrees of freedom. We will identify configurations 
related by the U ( n ) action 

Z^-UZU^ and p^-Up with U &U(n) (3-11) 

We further require that Z and p satisfy the matrix constraint 4 , 

irp [Z, Z^] + pp 1 = k l n (3-12) 

This constraint is the moment map for the action (3.11) with level k. We define the 
moduli space M n through the symplectic quotient, 

M n = {Z, p such that nfi[Z, Z^] + pp = fc} / U(n) 

This space has real dimension dim(A4 n ) = 2 n. The string theory construction of [42] 
shows that this space is related to the vortex moduli space 

Mn = Mn 

These spaces are conjectured to be isomorphic as complex manifolds, and have the 
same Kahler class. To our knowledge, there is no direct proof of this conjecture beyond 
the string theory construction provided in [42], 

4 As an aside: for relativistic vortices, the right-hand side of (3.12) is 27r/e 2 , where e 2 is the gauge 
coupling constant. Comparing the vortex equations (3.1) to their relativistic counterparts shows that 
this becomes k in the non-relativistic context. The fact that this is integer valued for vortices in the 
Chern-Simons theory will prove important below. 
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The matrix description provides a different parametrisation of the vortex moduli 
space. When the vortices are well separated, Z is approximately diagonal. The posi¬ 
tions of the vortices are described by these n diagonal elements. (The normalisation 
of 7 Tfi in (3.12) is associated to the vortex size.) However, as the vortices approach, 
Z is no longer approximately diagonal, reflecting the fact it is better to think of the 
locations of the vortices as fuzzy, spread out over a disc of radius (3.4). This feature is 
captured by the matrix description of the vortex moduli space. 

The moduli space M n inherits a natural metric through the quotient construction 
described above. This does not coincide with the metric on the vortex moduli space 
A4 n described in Appendix B. Nonetheless, there are now a number of examples in 
which computations of BPS quantities using A4 n coincide with those of computed 
from the vortex moduli space A4 n because they are insensitive to the details of the 
metric (see, for example, [43, 44, 45, 46]). Here we will ultimately be interested in 
holomorphic wavefunctions over the vortex moduli space. Assuming the conjectured 
equivalence of the spaces as complex manifolds, it will suffice to work with the matrix 
model description of the vortex moduli space. 

The Matrix Model Action 

It is now a simple matter to write the vortex dynamics in terms of these new fields. 
We introduce a U(n ) gauge field, ao, on the worldline of the vortices. In the absence 
of a harmonic trap, the low-energy vortex dynamics is governed by the U(n) gauged 
quantum mechanics, 

vortex j dt iir/iTv (Z^V 0 Z) + i(p^V 0 (p — k Tr oo (3.13) 

where T> 0 Z = doZ — i[ao, Z\ and V 0 (p = do — iao^p. The quantum mechanical Chern- 
Sirnons term ensures that Gauss’ law for the matrix model coincides with (3.12). This 
means that this action describes the same physics as (3.5). 

The action (3.13) is the quantum Hall matrix model , previously proposed as a descrip¬ 
tion of the fractional quantum Hall effect by Polychronakos [19] and further explored 
in [20, 21, 47, 48, 49, 50]. The connection to first order vortex dynamics was noted 
earlier in [23]. 

We note in passing that we’ve used the D-brane construction of [42] in a fairly indirect 
way to derive the quantum Hall matrix model. A more direct D-brane derivation of 
the matrix model was provided previously in [51]. It would be interesting to see how 
this work, or the string theory construction of [53], is related to the present set-up. 
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We would also like to add the harmonic trap to the matrix model. This too was 
explained in [19]. Spatial rotation within the matrix model acts as Z —» e ld Z , with the 
associated charge J = tt/iTy Z^Z. Adding this to the action, we get the matrix model 
generalisation of (3.8), 

'S vortex j dt m\i Tr (Z^VqZ) + i(p^V 0 (p — JcTy cio — loh/iTy (Z^ Z) (3.14) 

In the rest of this section, we describe the properties of this matrix model. Much of 
this is review of earlier work, in particular [19] and [20, 21]. However, we also make a 
number of new observations about the matrix model, most notably the computation of 
the charge and statistics of quasi-hole excitations. 

The Classical Ground State 

In the presence of the harmonic trap, the classical equation of motion for Z reads 

iD t Z — uZ (3.15) 

There is a unique time independent solution, with Z — 0, obeying [a 0 , Z) = ojZ. This 
can also be viewed as the statement that rotating the phase of Z is equivalent to a 
gauge transformation. There is a unique solution to this equation and the constraint 
(3.12) given by [19], 


with a 0 



(0 1 
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( ° \ 

Zo = ,/A 

0 y/2 



and (fi 0 = Vk 

0 

V 
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\V^J 

u diag(n — 1 

,n-2 ,..., 2, 

1 ,0). 





(3.16) 


As promised, Z$ is not approximately diagonal. This reflects the fact that individual 
vortices do not have well-defined positions. Nonetheless, we can reconstruct a number 
of simple properties of the vortex solution from this matrix. The radius-squared of the 
disc can be thought of as the maximum eigenvalue of ZqZ 0 [19]. To leading order in 
the vortex number n, this gives 



7r/i 
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which agrees with our the radius of the classical vortex solution (3.4). Meanwhile, the 
angular momentum of a given solution is J = Tr Z^Z. The angular momentum of the 
ground state is 

Jo = tt/z Tr (4Z 0 ) = knijl ~ 1} (3-17) 

which, to leading order in 1/n, agrees with the angular momentum of the classical 
vortex solution (3.9). 

The Quantum Ground State 

The quantisation of the matrix model (3.14) was initiated in [19] and explored in some 
detail in [20] and [21]. The individual components of the matrix Z and vector p are 
promoted to quantum operators, with commutation relations 

71" h [^abi ^ c d\ ^ad^bc and \jpat J &ab 

We choose the vacuum state |0) such that Z ab |0) = <^|0) = 0. However this does not, 
in general, correspond to the ground state of the theory because the physical Hilbert 
space must obey the quantum version of the Gauss’ law constraint (3.12). It is useful 
to view the trace and traceless part of this constraint separately. The trace constraint 
reads 

n n 

Twt = kn => VaVa = {k- 1 )n (3.18) 

a= 1 a— 1 

This means that physical states must have (k — 1 )n ^-excitations. Note that the 
ordering of the original constraint has resulted in a shift k —> k — 1. This will prove 
important below. 

Meanwhile, the traceless part of the constraint (3.12) tells us that physical states 
must be SU (■ n) C U (■ n) singlets. We can form such singlet operators out of Z^ and </d 
either from baryons or from traces. The baryonic operators are 

e a,..a„ ( ^t z t Pl)ai __ Pn)an 

where pi,.. ■ p n are, necessarily distinct, integers. The trace operators are 

Tr (Z tp ) 

There can be complicated relations between the baryonic and trace operators; explicit 
descriptions for low numbers of vortices were recently presented in [54], 
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The trace constraint (3.18) means that physical states contain exactly k — 1 baryonic 
operators. The harmonic trap endows these with an energy proportional to the number 
of Z' excitations, 


H = uj = W7T/J ^2 Z lb Z ba 
a,b= 1 

To minimise this energy, we must act with k — 1 baryonic operators, each with pi = i — 1. 
This results in the ground state 

| ground)*; = |0) (3.19) 


The angular momentum of this ground state coincides with that of the classical ground 
state (3.17). 


There is a close resemblance between these ground states and the Laughlin states [8] 
for n electrons at filling fraction v — 1/k, 


|Laughlin)& = \\{z a ~ z b ) k e E|z “ |2 

a<b 


'e ai - a "z°, z, 


a i 2 


. . . Z, 


n—11 


■EM 2 


(3.20) 


A formal map between the states was suggested in [20]. A more rigorous study, which 
we now review, was provided in [21] (see also [48]). The first step is to identify ap¬ 
propriate coordinates on the phase space of the matrix model. There are a number 
of different choices, none of which have preferred status. Here we use the coherent 
state representation suggested in [21]. We diagonalise Z and use its eigenvalues as 
coordinates on the phase space. (Non-diagonalisablc matrices have zero measure.) 

The first result of [21] is that the k — 1 ground state, |ground)i = |0), is precisely 
the v — 1 Laughlin state describing a filled Landau level, 


|0) = | Laughlin)! 


For k > 1, the map to the Laughlin wavefunction is not exact. Instead, the wavefunc- 
tions agree only at large separation 


|ground)fc —>■ |Laughlin)for \z a — z b \ ^ l/np 

However, the matrix model states |ground) ^ differ from the Laughlin states as the 
particle approach: the wavefunctions still vanish as z a —> z b , but not with the same 
power. 
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As we mentioned, there is nothing privileged about the choice of coordinates used 
above. One may try a different set of coordinates and see if there is better short- 
distance agreement with the Laughlin wavefunction. Indeed, other coordinates were 
suggested in [21, 48], although none of them provide an exact match to the Laughlin 
wavefunctions. 

The connection to vortices sheds some light on this. Because vortices are extended 
objects, there is no “correct” way to specify their positions as they approach. Corre¬ 
spondingly, it is not obvious that their physics is captured by a wavefunction describing 
point particles. Instead, the important questions are those which are independent of 
the choice of coordinates. The fact that the long-distance correlations in the matrix 
model ground states (3.19) coincide with those of the Laughlin wavefunction suggests 
that these states describe the same universality class of quantum Hall fluids. In the 
rest of this paper, we show that this is indeed correct. We show that excitations of the 
matrix model describe chiral edge modes and quasi-holes. In particular, the latter have 
charge 1/k and fractional statistics, in agreement with the excitations of the Laughlin 
wavefunction. 

3.4 Edge Modes 

The classical excitations of the matrix model were described in [19]. There are edge 
excitations of the droplet and there quasi-hole excitations although, for finite n, there 
is no clear distinction between these. There are no quasi-particle excitations which, 
given the spacetime picture in terms of vortices, is to be expected. We first study the 
edge modes and show that they form a chiral boson. 

The linear perturbations of the solution (3.16), consistent with the constraint (3.12), 
were given in [19]: they are remarkably simple, 

5iZ = (Zq)^ 1 and <5^ = 0 with l — (3-21) 

These were interpreted in [19] as area-preserving deformations of the disc, restricted to 
the first n Fourier modes. 

We now show that the dynamics is that of a chiral, relativistic boson. To do this, we 
write 


Z{t) = z 0 + YtCiWt 1 - 1 

1=1 
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with complex coefficients q. Plugging this ansatz into the action (3.14), we have the 
following expression for the effective dynamics of q, 


n r. 

s = vr/i I dt iT±(Zo~ lz o P ~ 1 ) c * d P + 

i , p = 1 ' 

where we have dropped the constant contribution (3.17). We need to compute two 
traces, both involving Z$ given in (3.16). The first is 

U- 1 

n/j Tr Zq~ 1 Zq~ 1 = QiSip with 0/ = —j—n(n — 1)... (n — l + 1) 

The second trace involves «o and can be readily computed by invoking the relationship 
l oZ 0 = [a 0 , Z 0 \, to give [a 0 , Zq P ] = —puZl P . The action for the perturbations can then 
be written in the simple form, 


Tr(a 0 [Z { 


tp-i 
o i 


rvl — 1 
Z 0 


]) - wiv(Zj-'Z*”- 1 ) 


etc, 


i 


s = E e ' 

i=i 


dt ( ic * Ci — l ole*c{) 


(3.22) 


This is the action for a real, chiral boson, defined on the edge of the Hall droplet. We 
parameterise the perimeter of the droplet by a G [0, 27lR) with R given by (3.4). The 
continuum excitations then take the form 


c(cr, t ) 



ci(t) 


with c_i = c* 


Then the action (3.22) becomes 


S = — J dtda d t cd a c+ ( y coR)d a cd a c 

This is the form of the action for a chiral boson proposed in [17], now truncated to 
the lowest n Fourier modes. The action describes modes propagating in one direction 
around the disc with velocity v = uR. A previous derivation of the chiral boson edge 
theory from the matrix model was given in [50], albeit in a model with a different 
potential. It is unclear to us how that derivation relates to the one above. 

Note that as n increases, the radius of the disc (3.4) scales as \fn, while the number 
of Fourier modes increases linearly with n. The density of modes therefore scales as 
l/\/n, suggesting the existence of a continuum d— 1 + 1 dimensional limit as n —> oo. 
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The Noncommutative Description Revisited 

The original motivation for the quantum Hall matrix model was to provide a finite 
n regularisation of Susskind’s non-commutative approach to quantum Hall fluids [52], 
Taking the n —> oo limit of the matrix model, one can effectively drop the field ip and 
the constraint (3.12) becomes 

i V 2 i - 2 Th . eB ext 

• x 1 = ! “T = ! — 

We interpret this as a non-commutative plane. Expanding the action around the state 
(3.16) gives rise to a Chern-Simons theory on this non-commutative plane, with fields 
multiplied using the Moyal product [52], The perspective offered here shows that 
this non-commutative theory provides a hydrodynamic description of the dynamics of 
n —> oo BPS vortices. 

There is no harmonic trap introduced in the non-commutative Chern-Simons descrip¬ 
tion. Because it arises from the expansion around (3.16), all perturbative excitations 
of the theory are edge modes of an infinitely large disc, now consigned to asympto- 
tia. However, these perturbation excitations are not the end of the story. There are 
many other non-perturbative bulk excitations. These correspond to separating vortices 
or, as we will see in the next section, creating a hole in the fluid of vortices. The 
non-commutative Chern-Simons theory is capturing these modes. 

However, we have already seen a different description of these modes from the per¬ 
spective of the d = 2 + 1 dimensional spacetime picture: they are the gapless, lowest 
Landau level of an interacting boson that we saw in Section 2.2. It appears that the 
Chern-Simons theory on the non-commutative plane is an alternative description of 
this lowest Landau level physics. 

3.5 Quasi-Holes 

Let us now return to a finite droplet of vortices. While the infinitesimal perturbations 
of the droplet describe edge modes, one can also consider finite deformations. Of 
course, if we make a large enough finite perturbation, then the droplet will eventually 
fragment into its component vortices. However, there are deformations for which the 
droplet retains its integrity, but with a hole carved out in the middle. These are the 
quasi-holes of the quantum Hall effect. 
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There is a simple classical solution describing a quasi-hole placed at the centre of the 
vortex [19]. It arises by integrating the n th Fourer mode, 


Z 


( 0 



+ Q ^ 

0 \/2 + q [ 


0 yjn - 1 + q 


\ ^qe inu}t 


0 ) 


(3.23) 


This obeys the constraint (3.12) and equation of motion (3.15) with a 0 = <ndiag(n — 
1, n — 2,..., 2,1, 0) and ip = ip Q . 


This solution should be thought of as a deficit of magnetic field in the middle of the 
Hall droplet [19] (see also [55]). In other words, it is a quasi-hole. Using the maximum 
and minimum eigenvalues of Z^Z as a proxy for the inner radius R\ and the outer 
radius R 2 of this annulus, we find 


R\ « 


kq 

7T/1 


and 


Ri~ 


k(n + q) 


TCp 


which is consistent with the magnetic flux quantisation (3.2) if B remains constant 
for R\ < r < i? 2 - We can subject this interpretation to a further test. The angular 
momentum of the matrix model solution is given by 


J = 7r/i Tr Z 


kn 2 

+ knq 


But we can also compute the angular momentum of an annular vortex by the same 
kind of calculation we used in (3.9). We find 


J 


k 

2 


-r 2 


dr 2irr 3 B 


' Ri 


kn 2 

~T 


+ knq 


confirming the solution (3.23) as a classical quasi-hole. 


There are, presumably, more complicated classical solutions, describing quasi-holes 
displaced from the origin, rotating with frequency to. Rather than searching for these 
classical solutions, we will instead describe their quantum counterparts. 
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Quantum Quasi-Holes 

We claim that the quantum state describing m quasi-holes, located at complex coordi¬ 
nates rji, i — 1,..., m, is 


oc JJ det(Z t -rjj) |ground ) k (3.24) 

i=i 

where we have allowed for a normalisation constant. 

Let us first motivate this ansatz. Multiplying by det(zU ~rj') is equivalent to taking 
one of the baryonic operators in the ground state (3.19) and replacing each occurrence 
of <p^Z' p by (p'Z' p (Z — r)y. Under the coherent state map of [21], where the eigenvalues 
of Z are used as coordinates, this gives 

\Vi,---,Vm) k ->• IJ(^o-^)|Laughlin) fc 

a 

which is indeed the Laughlin wavefunction for quasi-holes. 

As we vary the positions rji, the resulting states \r)i, ..., r/ m ) are not linearly indepen¬ 
dent. This reflects the fact that these holes are made from a finite number of underlying 
particles. Nonetheless, for [rji\ < R, with R = y/kn/nju, the size of the quantum Hall 
droplet (3.4), we expect the state to approximately describe m localised quasi-holes. 
This interpretation breaks down as the quasi-holes approach the edge of the droplet. 
Indeed, the states degenerate and become approximately the same for any value of 
|77i| R. We’ll see the consquences of this below. 

In the presence of a harmonic trap, the states (3.24) are not energy eigenstates 
unless rji = 0. Nonetheless, it is simple to check that the time-dependent states, 
\e lut 7]i, ..., a lujt Tj m ')ki in which the quasi-holes orbit the origin, solve the time-dependent 
Schrodinger equation. In what follows, we will compute the braiding of the time inde¬ 
pendent states (3.24). 

In the quantum Hall effect, the quasi-holes famously have fractional charge and 
fractional statistics. We now show this directly for the states (3.24). We follow the 
classic calculation of [18] in computing the Berry phase accumulated as quasi-holes move 
in closed paths. However, there is a technical difference that is worth highlighting. In 
the usual Laughlin wavefunction, the overlap integrals are too complicated to perform 
directly. Instead, one resorts to the plasma analogy [8]. This requires an assumption 
that a classical 2d plasma exhibits a screening phase. 
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A second route to computing the braiding of quasi-particles is provided by the link to 
conformal held theories [15], where it is conjectured to be equivalent to the monodromy 
of conformal blocks. The primary focus has been on the richer subject of non-Abelian 
quantum Hall states. Different approaches include [56] and [57, 58], the latter once 
again relying on a plasma analogy. See also [59] for an alternative approach to braiding. 

We will now show that the matrix model construction of the quasi-hole states (3.24) 
seems to avoid these issues and a direct attack on the problem bears fruit. We compute 
the Berry phase explicitly without need of a plasma analogy. 

Fractional Charge 

We start by computing the charge of the quasi-hole under the external gauge held. 
To do this, we consider a single excitation located at r/ = re ie . We then adiabatically 
transport the quasi-hole in a circle by sending 6 —> 9 + 2ir. If the quasi-hole has charge 
gqn then we expect that the wavefunction will pick up the Aharonov-Bohm phase 0 
proportional to the magnetic flux <f> enclosed in the orbit, 

27r 2 /yr 2 

0(r) = <f>g QH = irr 2 B ext q Qll = —-—g QH (3.25) 

where we’ve used the value of B ext = 27qu/e computed in (3.3), with e the charge of 
a single vortex. There is a more direct expression for 0, arising as the Berry phase 
associated to the adiabatic change of the wavefunction, 

0(r) = -i^ dO k (rj\-^\r)) k (3.26) 

Our task is to compute this phase. From this we extract ggn- 

To do this, it will help to introduce some new notation. We define the states \£li)k, 
with l — 0,..., n — 1 

l«i>* = [«“■.. ■ ■ Wu. 

>.“"4 1 (v , zt) ta ...qtzt”- I K]^ 2 |o> 

Each of these is an eigenstate of angular momentum, with J — J Q + 7r/i(n — /). We 
can expand the quasi-hole state (3.24) in this basis as 

n —1 

\v)k oc ^(-// t ) z |D z ) fc 

1=0 
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Figure 3: The Berry phase for a single quasi-hole in n = 1000 vortices with k = 3. The 
phase 0 (solid, red) and the expected phase for a particle of charge — e/k in the field B ext 
(gray, dashed) are both plotted. 

Because the \Qi)k have different angular momenta, they are orthogonal. We write their 
inner product as 


A(/, k'j Sip 


In terms of these inner products, the Berry phase (3.26) is simply written as 


@(r) = 27 ri 


EIU ***(*; fc)r M 

ELo Kl-,k) r* 


The computation of Ml: k ) is not straightforward. (Indeed, this is the step in the usual 
calculation where one resorts to the plasma analogy.) We find the following result: 


A (l-k) 



n—l—1 

n ( ka + x ) 

a =0 


k (ground | ground) k 


We relegate the proof of this statement to Appendix C. 


(3.27) 


Rather remarkably, the resulting sum can be written in closed form. We find 

n iFi(l — n, 2 — n — 1/k, nfiM/k) 

(n — 1 )k + 1 iFi( —n, 1 — n — 1/k, Tifir 2 /k) 

This is the ratio of confluent hypergeometric functions of the first kind. 
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The result (3.28) is plotted in Figure 3 for n = 1000 vortices and k = 3. The 
plot shows clearly that, for r < R , the Berry phase 0 coincides with the expected 
Aharonov-Bohm phase (3.25) if the charge of the quasi-hole is taken to be 

gQH = ~k 

This, of course, is the expected result [ 8 , 18]. 

Our Berry phase computation also reveals finite size effects. The magnitude of the 
Berry phase reaches a maximum of 2tt n at r = R , the edge of the droplet. Outside 
this disc, the Berry phase no longer increases and the picture in terms of quasi-holes 
breaks down. One can also use the result above to determine the size of the edge effects; 
numerical plots reveal them to be small as long as k <C n. 

There is another interpretation of the quasi-hole state (3.24): it is an excitation of 
the fundamental boson <fi in the Hall phase (2.17). Now the Aharonov-Bohm phase 
arises because this particle has charge 1 under the statistical gauge field with magnetic 
field B = —2'K/i/k. This is a pleasing, dual perspective. The vortices are solitons 
constructed from q k But, equally, we see that we can reconstruct 0 as a collective 
excitation of many vortices! 


Fractional Statistics 

Let us next consider the statistics of quasi-holes as they are braided. To do this, we 
consider a state with two excitations, 1771 , 772 )*,. It is simplest to place the first at the 
origin, 771 = 0, and transport the second in a full circle. This is equivalent to exchanging 
the quasi-holes twice and computes double the statistical phase. Of course, there is 
also a contribution from the Aharonov-Bohm phase O(r) described above and we must 
subtract this off. The resulting statistical phase, 0 stat , is then given by 

/■2tt o 

2 0 stat( r ) = J d6 _|o, rj)k - O(r) 
where, again 77 = re ld . 

To compute the statistical phase, we need yet more inner products. We define the 
states 


|SM* = [e" . ... (<?'z'XX'z" + 2 U» ■ • • (^z ,n+1 )«» 

e 11 . ‘•ftX , z%...(^z<’'-X, 
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Figure 4: The statistical phase for a quasi-hole encircling a second quasi-hole at the origin 
for n = 5000 and k = 3. The Berry phase 0 stat (solid, red) is plotted, together with the 
expected phase for a particle of statistics n/k (gray, dashed). 


This is similar to \£li)k, defined previously, except now each factor of Z' has been 
increased by 1. This is the effect of placing the extra quasi-hole at the origin. (For 
more general locations of the quasi-hole, we would need the obvious generalisations of 
these states |fij/,j)fc-) The states are again orthogonal. This time, we find the 

norm is given by 


k{^0,l\^0,l)k 

k (ground | ground) k 


(tt/x) 


l—2n 


n 


n—l—l 


a =0 


X 


7-1 


i I (ka + 1) 


El a +-*-) 

72—1 

(ka + 2) 


a =0 


a=l 


(3.29) 


With these functions, it is straightforward to determine an expression for the statistical 
phase in terms of a sum over n states. Once again, this sum has a closed form, this 
time given using regularised hypergeometric functions by 


20 stat/ r N = 2?r V 2 ( 2 -^ 2 (1 + l/fc,l-n;l + 2/fc,2-n-l/fc; 7r/xr 2 /fc) \ _ Q 

k y 2 ^ 2 ( 1 /^, —2//c, 1 — rz — 1/A:; J 

We plot this for n = 5000 and k — 3 in Figure 4. All other plots with k -C n have 
similar features. We see that there is clearly an intermediate, parametrically large 
regime, in which the pair of particles are both far from the edge of the disc and far 
from each other, where their exchange statistics are given by 

7T 

Qstat _ _ 
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This is the expected result for a quasi-hole at filling fraction v — 1 jk. 

4. Future Directions 

Supersymmetry has long proven a powerful tool to understand physics at strong cou¬ 
pling in relativistic systems, ft is clear that if this power could be transported to 
the non-relativistic realm, then supersymmetry may be employed to say something 
interesting about open problems in condensed matter physics. 

In this spirit, there have been a number of recent papers in which mirror symmetry 
(which can be viewed as an exact particle-vortex duality in d — 2 + 1 interacting 
systems) has been explored in the presence of external sources. This has been used to 
study impurities [60, 61], non-Fermi liquids [62] and the physics of the lowest Landau 
level [63]. It would be interesting to follow the fate of mirror pairs (or Seiberg duals) 
under the non-relativistic limit. 

The Laughlin physics described in this paper is, of course, just the tip of the quantum 
Hall iceberg. A long-standing open problem has been how to generalise the quantum 
Hall matrix model of [19] to more general filling fractions such as the Jain heirarchy. 
(See [64] for an attempt.) The perspective offered in this paper suggests a route. It 
is known that the most general Abelian quantum Hall state can be captured by the 
K-matrix approach [65], with an effective field theory given by several coupled Chern- 
Simons fields 

L = l-Kue^'AfaAl + + ... 

It is a simple matter to generalise this to a non-relativistic supersymmetric theory. 
However, the dynamics of vortices in these theories have not been well studied. A 
matrix model for the vortex dynamics in these theories would presumably furnish a 
description of the most general Abelian quantum Hall states. (A matrix model for 
vortices in a class of theories with product gauge groups was proposed in [61, 66].) 

Another obvious generalisation is to look at vortices in non-Abelian U(N) gauge 
theories. These were introduced in [42, 67]. The vortices now have an internal degree 
of freedom and the moduli space is given by 

N 

7T/i[Z, Z f ] + ^ = kl n 

i =1 

modulo U ( n ) gauge transformations. This model was previously discussed in the con¬ 
text of quantum Hall physics in [68, 69]. It would be interesting to better understand 
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what these states are and how they are related to the underlying non-Abelian Chern- 
Simons theories. 

Finally, we have restricted our attention in this paper to the case with chemical 
potential /i ^ 0. With fi — 0, the theory develops a non-relativistic super-conformal 
symmetry. An index was proposed in [70] but, beyond this, little seems to be known 
about the structure of these theories. With non-relativistic, conformal theories playing 
an important role in various cold atom systems, these super-conformal theories appear 
worthy of further study. 

We hope to return to some of these issues in the future. 
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A. Appendix: Non-Relativistic Limits 

The Lagrangian (2.1) can be derived by starting from a relativistic Chern-Simons mat¬ 
ter theory, with J\f = 2 supersymmetry, and taking a limit in which anti-particles 
decouple [7]. A number of other non-relativistic theories with different gauge groups, 
and more supersymmetry, have been constructed in this manner [71, 72, 73, 74], 

In this appendix, we review this non-relativistic limit. We construct a more general 
theory than that of [7], with gauge group U(N C ) and Nf matter multiplets transforming 
in the fundamental representation of the gauge group. We also show how the chemical 
potential term /j can arise in this limit. 

We restrict our attention to the bosonic fields and, only at the end, describe the 
generalisation to the fermions. The bosonic Lagrangian for the Af = 2 supersymmetric 
U(N C ) Yang-Mills Chern-Simons theory is 

Srei = - I d 3 x -j-^Tr (F, U F^) + Tr (A,d u A p - | A,A U A P ) + ^Tr (V,a) 2 

Nj 2 2 

+ \ V M + <t>W<f>i + y Tl '( ^ ~ l ’ 2 ) (A- 1 ) 

i =1 i 

Here a is the real, adjoint scalar which accompanies A M in the vector multiplet, while 
(f) t are fundamental scalars that live in chiral multiplets. We have included a Fayet- 
Iliopoulos term v 2 , but not real masses for the </>j. This can be done and results in 
different inertial masses in the non-relativistic limit. 

Before proceeding, it’s useful to perform some simple dimensional analysis. We work 
with h = 1. This, of course relates energy to inverse time scales. However, as we are 
ultimately interested in non-relativistic physics, we retain the speed of light c. This 
means that we have two dimensionful quantities, length L and time T. 

The factors of c in (A.l) are currently hidden in the notation. The measure is 

d 3 x = cdtd 2 x (A.2) 


while the derivatives are 

c z 

with a = 1,2 indexing spatial directions. Similarly, Aq = A t /c. The action is di¬ 
mensionless. The other fields have dimensions [A] = T _1 and [A a \ = [cr] = L~ l and 
[0] = L~ 1//2 . The parameters have dimension [k\ = 0 and [e 2 ] = [v 2 ] = L -1 . 
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We first take the infra-red limit, e 2 —> oo, to remove the Yang-Mills term. This also 
imposes the D-term as a constraint: 

^ = Ml - 1,2 ( A - 3 ) 

i 

Using this to integrate out the adjoint scalar a, the scalar potential terms in (A.l) 
become 


v =(t ) Z> M 2 ( A - 4 ) 

This kind of sextic potential is standard in supersymmetric Chern-Simons theories. The 
next step is to take the non-relativistic limit by discarding anti-particle excitations. To 
this end, we make the ansatz 


<j>i(x,t)= *— <j>i{x,t)e imcH (A.5) 

\J 2mc 

Here m is the mass of </>, which we read off from the quadratic term in the potential 
(A.4) 


2nv 2 


The key point of the non-relativistic limit is that <f> varies much more slowly that the 
frequencies me 2 set by the mass gap. In particular, this means that the ansatz (A.5) 
prohibits anti-particle excitations which scale as e +tmc2t . Plugging the ansatz (A.5) into 
the kinetic terms gives, after an integration by parts, 

\\T> t (j)\ 2 = — ( i|D,0| 2 + 2imftV t 4> + m 2 c 2 10| 2 ^) 

c z 2mc \c J 

The overall factor of 1/c is cancelled by the factor of c in the measure (A.2). The third 
term, m 2 c 2 \(j)\ 2 is designed to cancel the quadratic term in the potential. We now take 
the non-relativistic limit c —y oo. In doing so, we’re left only with the term linear in 
time derivatives. We can repeat this for all other terms in the action. In particular, 
taking a similar scaling of the potential (A.4), leaves us only with the quartic coupling 
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The same scaling can be applied to the fermions in the original Af = 2 theory. The end 
result is a U(N C ) Chern-Simons theory, coupled to Nf fundamental matter multiplets. 
To describe it, we revert to the notation (j> —> <f>. The final non-relativistic action is 

r u o i 

S= dtd 2 x + i'&l'Dtipi ~ —Tr e^ p (A fl d u A p - —ApA u A p ) 

i =1 

1 

y, {v a (t)\v a <j)i + ValpjValpi + 4BA) 

i —1 

[(tjfc )+ 2 (^tVj)(^Vi) ( A -6) 

For f/(l) with iVj = 1, this is the action (2.1) when /i = 0. (We have used the notation 
<9 0 rather than d t in (2.1).) 

The action (A.6) is invariant under superconformal transformations described in [7] 
and, in more detail, in [70]. The generators of the supersymmetries are the obvious 
non-Abelian generalisations of (2.9) and (2.10). 

Adding a Chemical Potential 

The action (2.1) also includes a chemical potential [i which plays a crucial role in our 
quantum Hall story. It is straightforward to add an analogous to term to the relativistic 
Lagrangian (A.l). It is 

Cp = /iTv (A 0 - a) (A.7) 

Obviously this breaks d = 2 + 1 Lorentz invariance. It preserves two of the four 
supercharges. Indeed, such terms are well known in the context of quantum mechanics 
models with A f = (0,2) supersymmetry and were first introduced in [75]. In taking 
the infra-red limit, the a term in (A.7) gets replaced by through the constraint 

(A.3). The resulting interaction terms of the non-relativistic theory are 

v = ^ 

Despite the fact that the relativistic theory with the deformation (A.7) preserves only 
one complex supercharge, both supercharges (2.9) and (2.10) are recovered after taking 
the non-relativistic limit. However, as we have seen, only Q 2 remains a symmetry of 
the spectrum. 
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B. Appendix: The Geometry of the Vortex Moduli Space 

In this appendix, we review a few basic facts about the geometry of the vortex moduli 
space. Suppose that we have at our disposal the most general solution to the vortex 
equation with winding n, 


<f>{x\X) and A z (x;X) 

We define 2 n zero modes ( 8 a qf>, 8 a A z ) to be the infinitesimal deformations which take us 
from one solution to another. 

8 a <t> = + ioi a (j) and 8 a A z = + d z a a (B.l) 

Here a a (x] X) is an accompanying gauge transformation. By construction, these zero 
modes solve the linearised versions of the vortex equations (3.1) for any choice of 
a(x,X). This ambiguity is fixed by further requiring that the zero modes obey the 
background gauge condition, 

2 ft 

d z 5 a A z + d~ 8 a A z = — (i(j)5 a (jJ - i^8 a ^) (B.2) 

The metric on the vortex moduli space M n is constructed by taking the overlap of the 
zero modes 

g a b = J d 2 x ^ (8 a A z 8 b A z + 8 a A z 8 b A z ) + (<5 a 0 t 5 b (j) + 8 a 4> 8 b ft) (B.3) 

In relativistic theories, this metric plays an important role: the low-energy dynamics 
of the vortices is described by a sigma-model on A4 n with metric g ab . The metric is 
known to be free of singularities. It is also Kahler, inheriting its complex structure 
from the natural action of complex conjugation on the fields. The associated Kahler 
form is 

^ab = i J d 2 x ^ (8 a A z 8 b A z - 8 a A z 8 h A z ) + (<S a 0 f 8 b <f> - 8 a (j) <y> f ) (B.4) 

We now show that this Kahler form governs the first order dynamics of vortices in our 
model. We will prove that the effective action for vortices is given by 


S , 


vortex 


dt X a {X)X a with dJ = O 


This result was previously derived in [24, 38, 39]. 
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We work in the usual spirit of the moduli space: we promote the collective coordinates 
of the static solutions to be time dependent: X a {t). We then substitute this time- 
dependent ansatz into the kinetic terms of the action (2.1). This results in an effective 
vortex action, 

S = J d 3 x ^ (A- Z A Z - A Z A,) + l - = f dtX a (X)X a (B.5) 


with 


FaiX) 



dAz 

OX 1 


■ A z ) + 




Note that the kinetic terms in (B.5) contain time derivatives rather than covariant time 
derivatives. This is because the A 0 terms in (2.1) multiply Gauss’ law and so necessar¬ 
ily vanish. Correspondingly, the expression for T„ above contains partial derivatives 
of holds which differ from the zero modes defined in (B.l) as they are missing the 
contribution from the gauge transformation a a (x]X). 

The 2-form f1 = dJ 7 is 

^ _ dT a dT b _ . f 2 k (dA z dA z dA z dA z \ | / d(fA d(p d(p dfl \ 

ab ~ dX b dX a ~ 1 ] X 7r \ dX a 0X b dX a dX b ) + \ dX a dX b dX° dX b ) 


Our goal is to show that Q a b = O a b , the Kahler form dehned in (B.4). The expressions 
look similar. They differ because the expression for fincludes extra contributions 
from the gauge fixing terms. We now show that these terms vanish. 

The proof is very similar to that given recently in [76] in the context of first order 
motion on the instanton moduli space. We take the difference 

~ /' 2 k f dA z dA z \ f. dcfA + dcj) \ , . 

n ab - n ab = I j d X - [g^d z a b - g^d- z a b j + J ~ (a++ b) 

= - I d 2 x a b ^+ <^0) - ( a +* b ) 

where we have integrated by parts to get to the second line. But the term in brackets 
vanishes, courtesy of the Gauss’ law constraint (3.1). We learn that dJ 7 = Q, the Kahler 
form, as advertised. Note that the proof above did not need us to use the background 
gauge fixing condition (B.2). While the metric (B.3) is sensitive to the background 
gauge condition, the Kahler form (B.4) is not. 


37 



C. Appendix: Overlap of Matrix Model States 


Our derivation of the fractional charge and statistics of quasi-holes relied on expressions 
for the norms of matrix model states given in (3.27) and (3.29). These results have 
been derived previously, most notably in the context of the Calogero-Sutherland-Moser 
model. Because these results are stated in a slightly different language, we use this 
Appendix to explain the connection. 

The quantum Hall matrix model is well known to be equivalent to the bosonic inte¬ 
grate Calogero-Sutherland-Moser model [19, 77]. This describes identical particles in 
one spatial dimension, placed in a harmonic trap and interacting via a specific inverse- 
square potential. To see the connection we begin, following [22], by working with a 
coherent state representation of all matrix model states. Firstly, expand Z = X\ +iX 2 
into Hermitian and anti-Hcrmitian parts, and let the overcomplete states \X , 0) be 
defined by 


Xi IA 0) = X \x, 0 ), 0 \x, 0 ) = 0 \X, 0) 


together with the normalisation 


\J<IX u \X.o) {.X ,, 


hj 


1 


where we have added hats to emphasise which symbols denote the quantum operators. 
With respect to these states, we can write all states in terms of a wavefunction by 
taking inner products with (X, 0|. This in turn gives us a way to compute the inner 
products of matrix model states by computing integrals over X, 0. In what follows, we 
work with convention 7T/i = 1. 


On these wavefunctions, 20 has the representation 


Z?- 2“ 1/2 

*7 


X ij ~ 


—) 

dx,J 


analogous to the raising operator of the more familiar Hcrmite polynomials. Hence, up 
to an overall normalisation, the states we are interested in all have wavefunctions of 
the form 


if(X, 4) = /(Zt) ... (^"-‘) a „] 


38 



where / is some homogeneous, gauge-invariant polynomial. Specifically, we have the 
following correspondence: 


|groimd) t : f(B) = 1 

l«M>* ■ f(B ) = det B- B^BX ■ ■ -B^U 

At given n, k, we will denote the state with a given choice of / simply by |/). 

One can evaluate the action of / on the state to obtain instead 

■MV <t>) = /( X) ... O?*”- 1 ) J*" 1 e-^ xa 

where at leading order f(B) ~ 2( deg /)/ 2 /(5). 

Now the relationship to the states of the Calogero model is seen by performing a 
change of variables: diagonalise X via X = UDU -1 , where D i3 = XiSij. Defining the 
Vandermonde determinant 

A-."—<•..<* 

i<j 

one sees that the wavefunction becomes 


= /(D) ■ A‘-‘ e-4* 2 . n.(/D)y 

Note that f(D) = f(x) is simply a polynomial in x, whose leading behaviour we can 
determine from /. Also, we can see that U, (j) have decoupled from x. 

Hence, taking account of the Jacobian A 2 for our change of variables, at a given n, k 
all inner products satisfy 

(f\g) = c n ,k J d n x e~ x2 A 2k f(x)g(x) (C.l) 

where c n ^ is a calculable constant which we do not need for our computation. 

As is shown in detail in [22], the key observation now is that the action of the matrix 
model Hamiltonian H on our wavefunctions is given by H = A -1 HcaiA, where iJcai is 
the Hamiltonian of the Calogero model at statistical parameter k. But the eigenstates 
of the Calogero model are known; they correspond precisely to the Hi-Jack polynomials, 
the multi-variable generalisations of the Hermite polynomials which are orthogonal with 
respect to the measure in (C.l). These are labelled by partitions A. 


39 



One may readily check that in fact f,g in equation (C.l) must be multiples of the 
generalised Hermite polynomials of Section 3 discussed in [78]. But now we can refer 
to Proposition 3.7 of that paper which is readily unpacked to give the ratios between 
the norms of general states. Concretely, their H\(x) have leading term 


H x (x) ~ 2' a ' 


and norms 

I H l(x) d(i(x) 


/dpi 


x 


( ... rp) 

\ x l '*'2 x r 


+ distinct permutations) 


number of distinct permutations 


_olA| TT ^a(c, d) + (a\(c, d) + l))(k(l\(c, d) + 1) + a\(c, d)) 

( ii A k(n-(c-l)) + (d-l) 


Here, |A| = JY Aj is the number of cells in the corresponding Young diagram, and 
a\(c,d),l\(c,d) are respectively the arm and leg length of the cell with coordinates 
(c, d) in that diagram. 


All that remains is to work out what choice of A and normalisation correspond to 
the examples of / given above for the matrix model states. It is easily found that 


| ground) : f = H ( „ 0 . 0) 

= f = .. 

/= 2 ( 2 ”- ,|/ 2 (”)ff( 2,2 .2,i,i.i) 

where there are n — l instances of 1 (resp. 2) in the second (resp. third) partition and 
then (3.27) and (3.29) both follow on evaluating the above product. 

It is very clear how this generalises to arbitrary states in the matrix model, especially 
if one realises the close relationship between the partition A and the original definition 
of the matrix model states \&i) k and |O 0 ,z) fc - 
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